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The phase coherence and pseudogap at the excitonic phase transition 
in a three-dimensional extended Falicov-Kimball model 
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Considering the phase coherence mechanism between the conduction band electrons and valence 
band holes we derive the excitonic phase stiffness for the three-dimensional strongly correlated 
electron system, ft is shown that the phase coherence between the conduction and valence bands 
is important for the excitonic condensation at the very low temperatures. We prove that the 
Bose-superfluid excitonic transition critical lines driven by the phase stiffness in the strong and 
intermediate coupling regimes differ considerably from the excitonic pair formation temperature 
transition lines. The true condensate transition temperature appears to be much smaller than the 
exciton pair formation transition temperature. We show that a BEC to BCS type crossover could 
be expected in our model. In the strong interaction case the system is driven into the insulator 
state with a BEC condensation at the very low temperatures, while in the small interaction region 
the condensate order is governed by the pair formation and in this case the system is evidently in 
the BCS-like regime. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



The Coulomb interaction between the conduction 
band electrons and the valence band holes causes in the 
solid state materials the formation of the new bound 
states of these two quasiparticles called the excitons. 1 
The low-density system of excitons behaves like a weakly 
nonidcal Bose-gas^ These new formations lead to the 
various interesting physical phenomena in solid state 
materials and they are the subjects of the intensive 
experimental and theoretical researches^— The ground 
state of the strongly interacting solid state system is 
unstable against the formation of the excitonic pairs 
and at the very low temperatures the new physical 
phases of the matter could develop. 1 Particularly, in the 
scenario of the semiconductor-metal phase transition a 
new phase develops approaching to the transition from 
the semiconductor side4 This state is called as the "ex- 
citonic insulator" (EI) state and is characterized by the 
strong coherence between the conduction band electrons 
and valence band holes. In the EI state a gap for the 
charged quasiparticle excitations is always present i It 
has been shown^ that the transition into the EI state is 
a second-order phase transition. Series of the recent ex- 
perimental investigations 3 -^ in TmSeo.45Teo.55 suggested 
the existence of the EI state in that material. Another 
example of the material with a well defined EI state is the 
quasi-one-dimensional Ta2NiSe5 with highly polarizable 
Se. The angle- resolved photoemission spectra on these 
compounds 1 -!, demonstrate that the ground state therein 
is an excitonic insulator. In general the solid state 
system in the EI state is considered as an electron-hole 
pair condensate 1 at sufficiently low temperatures. In 
the small interaction region or when approaching to the 



EI state from the metal part of the transition spectrum, 
the system is in the Bardeen-Cooper-Schrieffer— (BCS) 
superconducting condensate like regime with weakly 
bound electron-hole pairs and approaching from the 
semiconductor side the system shows typical Bosc- 
Einstein condensation 13 (BEC) behavior with the tightly 
bound excitons. Thus in the systems with the EI phase 
instability a BCS-BEC crossover— could be expected. 
This type of crossover mechanism is discussed in Ref. 
[TBI , where the authors find a transition from the BEC of 
excitons to a BCS state of electron-hole pairs. Using the 
screened Coulomb potential they derive the equations 
of quantum condensation using the real time Green 
function technique for a partially equilibrium system. 
In the context of the EI state, the three-dimensional 
spinless Falicov-Kimball model with the dispersive 
/-orbital electrons at half-filling has been analyzed 
recently*^ Particularly, from the behavior of the partial 
/ and c-electron densities of states, the authors find a 
crossover from a superconducting BCS-like condensate 
to a Bose-Einstein type excitonic condensate. A finite 
temperature transition critical line is established when 
the Hartree gap is non zero. The spontaneous symmetry 
breaking for the EI state and the BCS-BEC like crossover 
for the two-dimensional extended Falicov-Kimball model 
(EFKM) is discussed also in Ref. [It]- The authors study 
the EI state using the variational cluster approximation 
(VCA) and the Hartree-Fock (HF) approximation. It is 
shown that the coherence length of the excitonic pair 
spread wide for several hundred lattice spacing. This in- 
dicates that the excitonic insulator state typifies either a 
BCS-like weakly-bound condensate of electron-hole pairs 
or a Bose-Einstein condensate of preformed excitons. 
The nature of the electron-hole bound states is analyzed 
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in Ref. 18. Especially, at the semimetal-semiconductor 
transition the electron-hole pair condensation is found 
and the precursor effects for the electron-hole pair 
susceptibility is considered. 

At the temperatures below the temperature of the quan- 
tum degeneracy the excitonic quasiparticles obeying 
the Bose-statistics can spontaneously form the coherent 
states, which is the main characteristic of condensation 
in momentum space. At this level of treatment, to ex- 
plore the quantum physics of the excitons and in contrast 
with the mentioned theoretical works above ) 15-18 we 
derive the EI state by showing its direct relation to the 
phase coherence between the conduction band electrons 
and valence band holes. We show that this coherence is 
leading furthermore to the exciton Bose-condensation at 
the sufficiently low temperatures. Following the original 
papers , 19 ! 20 our model contains two type of fermions. 
The valence band /-electrons with the orbital energy 
e/ couple to the itinerant c-electrons with the orbital 
energy e c supposing that there is a mechanism for the 
f-c coherence. In the same spirit of the works ; 16 ' 17 we 
discuss the EFKM in the context of the path integral 
formalism and we derive the EI state solving the system 
of self-consistent set of equations. Basing on the form 
of the effective phase stiffness action derived in the 
model, we prove that the phase coherence between 
the c-orbital electrons and /-band holes is essential 
along the whole pair formation EI line. The important 
issue of the theory is the derivation of the excitonic 
condensate order. The temperature transition critical 
lines are derived employing the effective quantum rotor 
model. We show that temperature transition critical 
lines obtained for the Bose-superfluid condensation of 
excitons in the three-dimensional system of excitons are 
considerably different from the excitonic pair formation 
temperature transition lines. Hence, we show that 
pair formation transition and the excitonic condensate 
transition are two different lines which are in coincidence 
only in the very small Coulomb interaction region, where 
the system is driven to the BCS-like state. Our theory 
provides an alternative mechanism for the BEC-BCS 
crossover discussed earlier. 

The present article is organized as follows. In the Section 
[IT] we provide the Hamiltonian of the model EFKM, then 
in the Section [En] we introduce the new decoupling po- 
tentials and we handle with the four fermion interaction 
term in the Hamiltonian. In Section IIVI we derive the 
excitonic pairing gap parameter, single-particle excita- 
tion energies and the average particle concentrations for 
each band. In Section [V] we derive the phase stiffness 
parameter and the self-consistent equation for the 
excitonic condensate transition temperature in EFKM 
and we find the complete quantum critical solution of 
the EFKM model. 



II. THE HAMILTONIAN 

We consider the Hamiltonian of the extended spinless 
Falicov-Kimball model in which the electrons with one 
spin direction are localized, while those with the other 
spin direction are moving 

n = -tJ2 [c(r)c(r') + h.c] - (// - e c ) ^ n c (r) 

<rr'> r 

-tj2 [mm+h.c.} _( M _ e/ )£ n/ (r) 

(rr') r 

+tf5> c (r)n,(r), (1) 

r 

where (rr'} runs over the nearest neighbor sites, c(r) 
(c(r)) are the creation (annihilation) operators of the 
electron of the c-orbitals at the site with the position r, 
t is the hopping integral for the c-electrons and e c is the 
on-site energy level. Similarly, /(r) (/(r)) are the cre- 
ation (annihilation) operators of the /-orbital electrons 
and t is the hopping integral for the /-electrons. Further- 
more, e/ is the on-site energy level of the /-orbital. The 
chemical potential /i entering in Eq. (UJ is defined for our 
grand canonical system. We suppose that the chemical 
potentials of both orbitals are the same as in the work 
in Ref. \VA Parameter U, entering in the last term of the 
Hamiltonian, is the Coulomb repulsion between two types 
of electrons. n c (r) and n/(r) are the c- and /-electron 
density operators and they are defined as usual by the 
relation n x (r) — x(r)x(r) with x — c, / for the c and 
/-orbital electrons respectively. We will use throughout 
the paper the units ft = fcs = 1. We will use t = 1 as 
the measure of the energy. We consider also the follow- 
ing values for the band parameters e c — and ef = — 1 
also in accordance with the Ref. With this consid- 
eration the c- electrons are itinerant and the /-electrons 
are quasilocalized on the atomic sites. The equilibrium 
value of the chemical potential will be determined from 
the half-filling condition (the grand canonical approach) 
of the average particle density functions, i.e. n c = 1 — n/, 
where n x = (n x (r)) is the average particle density. 



III. DECOUPLING PROCEDURE AND 
GAUGE FIELDS 

We will decouple the four fermionic term in the Hamil- 
tonian Eq.([T|) within a path integral formalism. For this, 
we rewrite first the density product in the last term in 
Eq.([T]) in the equivalent form 



n c (r)n/(r) 



n 2 (r) n 2 (r) 



2 2 ' 

where we introduced the short-hand notations 

n(r) = n c (r) + n/(r), 
n(r) = n c (r) - n/(r). 



(2) 



(3) 
(4) 
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With the new notations we can rewrite the Hamiltonian 
Eq. (JlJ of the system as 

U = -t £ [5(r)c(r') + h.c] - p £ n(r) 

(rr'> r 

-t [mm + h.c] + ^ 

(rr'> r 

+uY,\b 2 ^)-^^)\- ( 5 ) 

r 

We put here p = ji — e and e = (e c + ey) /2 is the average 
band energy parameter. 



A. Functional integral formulation 

Dealing with the fermions, it is convenient to introduce 
the Grassmann variables c(rr) and /(rr) at each site r 
and imaginary time r varying in the interval < r < /3, 
where /3 = l/k B T (with T being the thermodynamic 
temperature). The variables c(rr) and /(rr) satisfy the 
anti-periodic boundary conditions x(rr) = — x(rr + j3). 
The partition function of the system of fermions written 
as a functional integral over the Grassmann field is 



Z = / [VcDc] / [DfVf] e 



S{c,cJ,f] 



where the action in the exponential is given in the path 
integral formulation as 

rP 

S[c,cJ,f}=S B [c,c}+S B [f,f} + drUir). (7) 

Jo 

Here S B [c, c] and Ss[f, f] arc the fcrmionic Berry terms 
for the c- and /-electrons respectively. They are defined 

as 



S B [x 1 x]^'^2J cLtx{yt)—x{vt). 



(8) 



Then we decouple the quadratic density terms in Eq.([5]) 
using the Hubbard-Stratanovich transformation and by 
introducing the new decoupling variables V(rr) and 
q(yt) conjugated to the electron density operators n(rr) 
and n(rr). For the quadratic term in the exponential of 
the partition function proportional to n 2 (rr) we have 



ex P (^-^J2j Q d ™ 2 ( rr )^) 



= / [W\ e 



(9) 



The decoupling of the quadratic term proportional to 
n 2 (rr) in the exponential of the partition function in 



Eq.© is also straightforward. We have 
rP it 



° xp (?r 



[Dq], 



dr—n (rr) j 
■E t /„ S dr\^p-- e {rr)h{vr)\ 



(10) 



Shifting the integration variables in Eq. (|10j) we will have 

s Er/>{-^-«W[«W-V]} 

(11) 



[Vq], 



The saddle-point evaluation of the obtained integral gives 
for g 



go = Un — 



(12) 



where n — (fi(rr)). As the result of the decoupling we 
will have the following "Zeeman-like" contribution to the 
partition function 



exp J dTn(rT)fi^j 



(13) 



(6) with the effective chemical potential fi e 



Un. 



After combining the exponential in Eq. §§§ with the effec- 
tive p chemical potential term in Eq.© we can decom- 
pose the variables V(rr) into a static and the periodic 
parts as 



V(rr) = V (r)+V(rr). 



(14) 



The integration over V(rr) variables now becomes the 
integration over the scalar static variables Vo (r) and the 
integration over the periodic fields V(r-r) 



[DV] ... = / [Wo] 



(15) 



For the periodic part in Eq. ([i"4")) we define the emergent 
Faraday-type relation introducing the U(l) topological 
phase field (^(rr) 



V(rr) 



dip(rT) 
dr 



cp(rr) 



(16) 



Thus, for the dynamic part, we transform the integration 
over the gauge variables V(rr) into the integration over 
the generic phase variables y(rr) in accordance with the 
relation in Eq.([T6|) 



VY 



— > 



[Dip] 



(17) 



In performing the integration over the phase field one 
should take into account that the field configurations sat- 
isfy the condition 



<p(rf3) — <^(r0) = 27r?7j(r), 



(18) 
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where m(r) = 0, ±1, ±2, ... are the U(l) topological wind- 
ing numbers, which are defining a set of the equivalent 
topological sectors. Thus the variables f{vr) only with 
discrete configurations are matter. The proper integra- 
tion measure in Eo.(|17p over ip variables is a multiply- 
connected topological domain and 



J r i m JO 

[m(r)] r 

rf(r/3)— (^(r)+27rm(r) 



rr 



(19) 



<p(rO)=<p(r) 



where in each topological sector the integration goes over 
the gauge potentials with the winding number equal to 
m(r). The periodic part, on the other hand, can be de- 
composed into the Fourier series as 



1 x r 



V(rr) = -^[v(ro; n )e i 



+ c.c. 



(20) 



with uj n — 2nn/f3 (n — 0, ±1,±2, ...) being the Bose- 
Matsubara frequencies. The decomposition of the charge 
field conforms with the basic m — topological sector, 
because 



d<p(rr) 



l * Or 



drV(rr) = 0. 



(21) 



For the scalar static part Vo(r) we have the following 
functional integral 



[dV ] e 



{-^+iVoW[n(rT)-I/p]| 



(22) 



The saddle-point value for the static part Vo(r) is given 
as Vq P = iUn — ip,, where n is the total average particle 
density n = n c + rif. And we have the contribution to 
the partition function in Eq.([B]) in the form 



[Dtp] 



,S[<p] 



(23) 



where the effective phase only action S[<p] is given as 



S[tp] =Y dT 



(p 2 (rT) n . 



- — (p{rr) — iip(rT)n(jT) 



(24) 



2U iU 



Combining the exponentials in Eq. (fT3)) and in Eq.([22l 
(after saddle-point evaluations) we will have our decou- 
pling result in the form 



cxp (^^2 J dT & n (. rT ) + Me^( rT )]^ J [VP] 



(25) 



where the effective chemical potential jx is defined as 
fi = Un — /2. The partition function of the system af- 
ter decoupling procedures will be 

Z = J [VcDc] [VfVf] [Dip] e ~s[c,cjj, v }^ ( 26 ) 

where the action S[c, c,f,f, ip] in the exponential is 
S[c,c, f, f,ip]=S[ip]+S B [c, c] + S B [f, f] 



-ty] I dr [c(rr)c(rV) + h.c] 

i i\J0 
(rr') 

-tJ2 / dr[f(rT)f(r'T)+h.c] 

+ 22 dT [/* n ( rr ) + M e n(rr)] . (27) 

r J 

We see that a strongly fluctuating imaginary term ap- 
pear in the expression of the phase action in Eg. (124)) ). 
In order to eliminate this term we will apply the U(l) 
gauge transformations to the fermionic variables c(rr) 
and /(rr). 

B. The U(l) gauge transformation 

We perform the local gauge transformation to the new 
fermionic Grassmann variables a(rr) and b(rr) in or- 
der to eliminate the last imaginary term appearing in 
the expression of the phase action in Eg . (|2~4"|) . Doing so, 
we cast the strongly correlated fermionic system into a 
weakly interacting one, which is submerged into the bath 
of strongly fluctuating gauge potentials (whose interac- 
tions are governed by a high energy scales of U) . For the 
c-orbital electrons the U(l) transformation is 



c(rr) 




c(rr) 





e Mrr) o 
g-^rr) 



a(rr) 
a(rr) 



(28) 



" /(rr) " 











' 6(rr) 


. /(rr) _ 







e ~i<p(rr) 




6(rr) 



where we introduced the new bosonic phase variables ip. 
For the /-orbital electrons the similar transformation is 



(29) 



After those transformations the electronic action of the 
system in Eg. (|27f takes the following form 

S[a,a, b, b, tp] = S [ip] + S B [a, a] + S B \b, b] 

rP 



(rr') 
(rr') 



dT 



dr 



a(rr)a(r , T)e- l [ v(rr) - ¥ ' (r ' T) ] + h.c. 



6(rT)6(r'T)e- l M rT) - v(r ' r) ] 



h.c. 



rP 

V / dr \fin(rr) + ll b u{vt)] 
^ Jo 



(30) 
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The phase action Sq [<p] in Eq. (f30|) is derived as 



Then action of the system takes the form 



s o[v] = X) / dT 



/3 



<^ 2 (rr) . 



-77^( rr ) 



(31) 



2U iU' 

The partition function of the system in new variables is 



Z = J [DaDa] [VbVb] [Dip] e 



— S[a,a,b,b,tp] 



(32) 



This form of the partition function is very general and 
furthermore we will use it for deriving the effective ac- 
tions for the fermions and for the phase sector. 



IV. RESULTS OF CALCULATIONS 

A. Effective fermionic action and propagators 

We start with the phase averaged expression of the ac- 
tion in Eq.© after the U(l) transformations given by 
the Eqs.(2J) and (25). Applying the Hartree-Fock (HF) 
like approximation (see Ref. 1231 ) to the four-fermion term 
in the Hamiltonian of Eq. (TJ , the phase averaged expres- 
sion of the total action can be derived. We present here 
the Fourier space evaluations (the derivation in the real 
space is very similar) . The Fourier transformation of the 
fermionic variables a(rr) and 6(rr) is given by 



(33) 



with x = a,b for the a and b type electrons. N is the total 
number of k-points in the first Brillouin zone (FBZ) and 
v n = n(2n + l)//3 are the Fermi-Matsubara frequencies 
with n = 0, ±1, ±2, ... . We introduce the local excitonic 
gap parameter A as in the Ref. [171 



A = U (a(rr)6(rr)) 



(34) 



Then, the HF decoupling of the interaction term in the 
Hamiltonian in Eq.(T]) is 

n a (rT)n b (YT) w (ra (rr)) n b (rr) + (n fe (rr)) n Q (rr) 
-iA(0, 0)a(rr)b(rr) - I A(0, 0)6(rr)a(rr), (35) 

where n a (rr) and nb(rr) are the electron density oper- 
ators after the U(l) gauge transformation. Due to the 
symmetry of the U(l) transformation they are practi- 
cally the same operators as n c (rr) and n/(rr), because 
of the relations n a (rr) = n c (rr) and n b (rr) = n/(rr). 



S [a, a, b, b] = ^ ak(v n )b k (v n ) 



777 £ O k (^„) (^cff - ilJ n - *k) ak(fn) 



ki/„ 



"77 £ &k(^n) (/4ff ~ - *k) &k(^n), (36) 



where A = A(0, 0) is the local excitonic gap parame- 
ter and the effective chemical potentials (j%g and n b cS are 
defined as 



Mcff = £ a - M + Un b + i (<p(rr)) , 
Meff = e fc - [i + Un a + i (<p(rr)) . 



(37) 
(38) 



n a and n& in the definitions of Eq. (l37l) and Eq.([38) are 
the average values for the respective density operators 
and n x = (n x (rr)) with x = a, 6. The renormalized 
hopping amplitudes tk and tk are given by the relations 



tk = 2i.ge (k) , 
4 = 2i> (k) , 



(39) 
(40) 



where 



e (k) = cos(d x k x ) + cos(d y k y ) + cos(d z k z ) (41) 

is the three-dimensional lattice dispersion relation with 
di, i — x, y, z, being the components of the lattice spacing 
vector d = r — r' with r and r' nearest neighbors sites. 
The bandwidth renormalization factor g is given as 



g = / e - i ['p(r-r)-ip(r'T)] 



|r-r'|=d 



(42) 



Introducing the vector space notations, we can rewrite 
the action in Eq.(35) in the form 

S [a,a,b,b] = — ^ ^ [ak(v n ),bk(y n )] S _1 (k^„) x 



bk(v n ) 



(43) 



Here Q 1 (kv n ) is the inverse of the Green function matrix 
given as 



a _1 (kz/„) = 



«kK) -a 
-A ££(!/„) 



(44) 



where single particle quasienergies £^(v n ) and ££(^„) are 
given after Eq. ([36|) as 



(45) 
(46) 
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The explicit expressions for the normal and anomalous 
excitonic propagators can be derived from the effective 
action in Eq. (j43|) . The normal propagator is defined as 
usual by (see Ref. HI) 



g xx (YT,v'T') = - (i(rr)i(r'r')) 



(47) 



The anomalous or the excitonic propagator is defined in 
the theory as 



£ ab (rT,rV) = (6(rr)&(rV)> 



(48) 



Then, using the standard functional derivation 

technique^ we have for the normal electron prop- 
agators 

i „ i [ k ( r - r ')- !/ »( r - r ')] 



k//„ 



fkK)^K)-|A| 2 ' 



i p i[k(r-r')- 1/ n( T - T ')] 



(49) 



The anomalous propagators are derived as 



Q ah { TT i T>T>) = —=-Y — 

\ ' ^ RM pa 



A e - i [ k ( r_r ') _1/ "( r_r ') 



£ bQ (rr, rV) 



A e - i [ k ( r ~ r ')-M r -' r ')] 

^v^^K)^K)-|Ap- 



(50) 



The local expressions of the normal and excitonic prop- 
agators will be used to find the system of self-consistent 
equations of the theory. The anomalous propagators of 
the form given in Eg. (1501) will determine the phase stiff- 
ness parameter (see Appendix iBl) . 



B. The excitonic gap parameter 

The general definition of the excitonic order parameter 
is terms of the initial c-band electron and /-band hole 
operators c(rr) and /(rr) is 

* = U lim (c(rr)/(r'r)) . (51) 

r— >r' 

Or after the U(l) transformations we have 

* = U lim (a(rr)&(rV)) / e *M^)-^( rV )] \ . (52) 



Thus we represented the excitonic order parameter fasa 
product of the U(l) gap parameter U (a(rT)b(r' t')) with 
a phase coherence transition probability function 



* = lim A(rr, r'r)Q 9 (r - r'), 



(53) 



where we introduced the excitonic gap parameter 
A(rr, r'r) as 



A(rr,rV) = U (a(rr)6(r'r)) 



(54) 



and the probability function G v {v — r') is, in fact, a phase- 
phase correlation function defined as 



Mr-r') 



(55) 



As we will see in the Section |Vj only a nonzero value of 
the gap parameter is not sufficient for the possible pre- 
diction of the excitonic condensation. To achieve this, 
also the phase variables must become stiff and coherent, 
which implies that Q v function defined in Eq. (|55[) is not 
vanishing at least for the fundamental state of the sys- 
tem in which the Bose particles are condensed in the 
k = momentum state of the reciprocal lattice space. 
This is the quantum state of the system, where a macro- 
scopic number of particles are situated in the same state 
of the reciprocal space. The Fourier transformation of 
the phase coherence factor G v (r — r') could be written 
as a sum of two components, one of which is describ- 
ing the condensate transition probability and the other 
one is just a usual fc-summation over the reciprocal wave 
vectors, defining the non-condensed part. Hence we can 
write 



?,(r - r') = <A 2 + ^ $>* k (- 

k/0 



where 



is the condensate transition probability with 



— / Jv{ rT ) 



(56) 



(57) 



(58) 



The non-zero value of the average in Eq. (|55|) is necessary 
for the transition to the Bose-superfluid state. The fact 
of the appearance of this non- vanishing expectation value 
leads to a phase transition from a symmetric phase tpo = 
to an unsymmetric phase ipo ^ 0. We can write now 
the excitonic order parameter as 



* = A 



^E^(k) 



k#0 



(59) 



where A is the local excitonic gap parameter A = 
U (a(rr)6(rr)) = A(0,0). * order parameter is consist- 
ing of the condensate order plus the on-condensate quasi- 
particle excitations. Indeed, the second term in Eq.([59|) 
gives contribution to the excitonic pairing mechanism 
without the condensate order. At the very low tempera- 
tures, at which the excitonic condensation could happens, 
we can neglect the on-condensate part in Eq. (|59j) and we 
have only the condensate order 



(60) 
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Indeed, the generic U(l) gauge field variables (p intro- lattice 
duced in Eq. p^j) . being apart the phase flux tubes at- 
tached to the fermionic variables, are also the main mech- 
anism for the Bose-supcrfluid state in the system. The 
exact prove of it is given in the Section [V] where the 
excitonic condensation is demonstrated using a quantum 
rotor model analogue to the 3D spherical model. 21 



C. Pairing gap 

We will show below how the excitonic gap parameter A 
could be calculated by solving a system of self-consistent 
equations for our system and, later on, in the Section [V] 
we will derive the condensate order and show how the 
phase stiffness mechanism is responsible for the temper- 
ature transition critical lines to the condensate order in 
the system. From the expression of the total action of the 
system in Eq. (|43p . we can derive the equations for the av- 
erage electron densities n a and rib corresponding to the 
a- and 6-orbitals respectively and also a self-consistent 
equation for the excitonic gap parameter A. Then we 
will have 



n a = £ aa (0,0), 
n b = g bb (O,0), 

A = UQ ab (0,0). (61) 

Using the analytical expressions for the propagators given 
in Eq. (|49l) and Eq.([50|). we can rewrite the system in 
Eq. (l61l) in an equivalent explicit form. We have 



jfY t [nF(E+) + n F (E^)]=l, (62) 
n F (E+)-n F (E k ) 



N 



£6. 



UA^n F {E+)-n F (E^) 



(63) 



(64) 



where wc used the fact that at half-filling n = n a + rib = 1 
and we defined n = n a — ri\,. n F is the Fermi-Dirac 
distribution function n F (e) = l/(e^ e + 1) and the energy 
parameters E£ and E^ are given by 



=2 [-^±\/S 



2 +4A 2 



(65) 



where Ak — t\c — /zjjg + tk — /i" ff and the quasiparticle dis- 
persion is £k = £k + Meff — ^ k — /^cff ■ The quantities n a , n&, 
A and [i can be determined by solving numerically the 
above equations in a self-consistent way. From the sum- 
mations over the wave vectors in Eqs. (|6"2"l) . (|6"3"|) and (|6~4")l 
we can pass to the integrations by introducing the appro- 
priate density of states (DOS) for the three-dimensional 



P3D 



(0 = ^£* 



C-^cos(kd) 



(66) 



For the simple cubic lattice geometry the density of states 
is given as 

1 ,min(l,2-C) 9 fl - IQ) 

P3D(0 = — I dy- ' ' 




(67) 

where O(x) is the Heaviside step function and K(x) is 
the elliptic function of the first kindj^The plot of the 
function psd{x) is presented in the inset of Fig.[TJ where 
the electronic density difference function n is plotted as a 
function of the normalized Coulomb interaction parame- 
ter U for the excitonic pair formation region, i.e. when 
A 7^ 0. Two different values of the hopping parameter 
t are considered and the solutions at T = K are pre- 
sented. It is clear from Fig. [1] that in the strong coupling 
limit, i.e. for the higher values of the Coulomb interac- 
tion, the system is in the insulator regime, because at 
the upper boundary of the Coulomb interaction parame- 
ter the 6-band is fully occupied (n.& = 1) and the a-band 
is totally free (n a = 0). It is also clear from Fig. [T]that 
with increasing the hopping integral the system reaches 
the "excitonic insulator" phase at the higher values of 
the Coulomb interaction. This is due to the fact that the 
higher hopping amplitudes are destabilizing the insulator 
phase and the limit n a = is achieved in the strong cou- 
pling case. In Fig. [2] the excitonic pairing gap parameter 
is plotted as a function of the Coulomb interaction energy 
in units of t and for different values of the hopping am- 
plitude t equal to t = —0.3 for the upper curve and equal 
to t = —0.1 for the lower curve. Decreasing (in absolute 
value) the hopping integral t, we observe that the corre- 
lation region becomes narrow, vanishing for the smaller 
values of U/t. The correlation amplitude also has been 
reduced. It is clear that in the three-dimensional case 
we have the long-range correlations between the electrons 
and holes. The excitonic gap is not vanishing for a rather 
large domain of the Coulomb interaction in difference 
with the results for the two-dimensional square lattice ge- 
ometry given in Ref. The authors in Ref. [T3, analyz- 
ing the EFKM Hamiltonian using the Hartree-Fock ap- 
proximation, derive the gap parameter for a narrow range 
of the Coulomb interaction strength U. The obtained 
values for the lower and upper bounds of the Coulomb 
interaction in Ref. [l7] are about (U c i, U C 2) — (0.66, 6.95), 
which differ considerably from our results especially for 
the large hoppings (see the upper curve on Fig. [5] for 
t = —0.3). In our case, the values for the lower and upper 
bounds of the normalized Coulomb interaction are given 
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2 4 6 8 10 



u/t 

FIG. 1: The average particle density difference function h as 
a function of the dimensionless interaction parameter U /t for 
two different values of the hopping integral t and at T = K. 
In the inset of the figure the density of states (DOS) for the 
three-dimensional cubic lattice is presented. 



2.5 




2 4 6 8 10 



u/t 

FIG. 2: The excitonic gap parameter A as a function of the 
normalized Coulomb interaction energy U (in units of t) at 
T — K for two different hopping parameters. 



with the intervals (U cl ,U c2 ) = (0.7,10.7) for t = -0.3 
and (U cl , U c2 ) = (0.8, 7.8) for the case t = -0.1. At the 
upper bounds of the Coulomb interaction we are in the in- 
sulator regime. In Fig. [3] the average density difference h 
is presented along the excitonic pairing boundary, where 
the pair formation begins, thus at the transition from the 
symmetric A = to the unsymmetric phase A ^ 0. We 
see that in this case we are starting nearly with the same 
values for the average densities n a and as in the case 
of the pair formation i.e. n a w 0.39 and rif, ks 0.61 at 
U = 0. They are different for the intermediate values 
of the Coulomb interaction parameter. It is clear from 
Fig-Elthat, approaching to the boundary at which A = 0, 
we are reaching the insulator state more drastically than 
in the case with the excitonic pairs. Indeed, for exam- 
ple at U/t — 6 we have n a = 0.107 and n& = 0.893 at 
the boundary A = 0, while, in the case when A / 0, 
we have for the same value of the Coulomb interaction 
n a = 0.1625 and rib — 0.8375. Thus, approaching to the 
excitonic pairing boundary, the average density of the 
valence band electrons is slightly increasing. The lines 
in Fig. [4] represent the corresponding chemical potentials 
/xa, along the excitonic transition boundary. The upper 
critical line corresponds to the case t = —0.3 and for the 
lower curve t — —0.1. The single particle excitation en- 
ergies at the pair formation boundaries are sensible to 
the changes of the hopping amplitude t. Smaller val- 
ues of t correspond to the smaller excitation energies. 
The chemical potential, in the case of the presence of 



excitonic pairs (derived from the system of equations in 
Eq. (l6lj) ). forms a well defined band and, consequently, 
implies the formation of the single particle excitation gap 
A g = fi max - pt mm , where fi max and [i mm are the upper 
and lower bounds of the chemical potential. The lower 
bounds of the chemical potentials for different values of 
the hopping integral t are presented on the Fig. [5] We see 
from Fig. [5] that the evolution of the chemical potential 
is continuous from weak to strong coupling limit. It is 
worth to mention that the similar behavior for the chem- 
ical potential is found in Ref. H^, where the authors ex- 
amine the Bose-condensation in an attractive interacting 
fermion gaz. In the inset of the Fig. [5] the average density 
difference n is presented as a function of the excitonic gap 
parameter normalized to the hopping amplitude. For 
the small values of the Coulomb repulsion the single par- 
ticle gap is totally degenerated and we have A s = 0. In 
this region the chemical potential is coinciding with the 
values presented on the Fig. 2) In a weak to intermedi- 
ate coupling regime, the single particle gap A g and the 
excitonic gap A both are increasing while in the strong 
couping limit (U > 8) the excitonic pairing gap parame- 
ter decreases rapidly with increasing U (see Fig. [5]) and 
the single particle excitation gap remains open. 



9 




FIG. 3: The average particle density difference function n 
along the pair formation transition boundary (A = 0) as a 
function of the dimensionless interaction parameter U/t. Two 
different values of the hopping integral t are considered. 
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-1 



- f=-0.3 

- f=-0.1 



■■■>■' 



10 



U/t 



FIG. 4: The chemical potentials along the pair formation 
transition boundary (A = 0) as a function of the dimension- 
less interaction parameter U/t. Two different values of the 
hopping integral t are considered. The chemical potentials 
are completely degenerated at the transition lines. 



FIG. 5: The chemical potentials accompanying the excitonic 
pair formation as a function of the dimensionless interaction 
parameter U /t and at T = K. Two different values of the 
hopping integral t are considered. Degeneracy of the chem- 
ical potentials is lifted due to the EI state. Only the lower 
bounds of the chemical potentials are presented. In the in- 
set, the variation of the density difference between two bands 
is presented as a function of the excitonic gap parameter for 
£=-0.3. 



V. PHASE COHERENCE TRANSITION LINES 

A. Phase stiffness at T = K 

In this section we examine the phase stiffness mecha- 
nism between the conduction band electrons and valence 
band holes. We will show how the phase coherence be- 
tween these two solid state bands is important for the 
excitonic condensation order at the very low tempera- 
tures. For considering the phase stiffness we are interest- 
ing in purely phase action and we will integrate out the 
fermions in Eq. ([521) . Then we obtain the effective phase 
action. The partition function after integration is 



where 



S cS [if] = - In / [VaDa] [VbVb] e 



(68) 



(69) 



Basing on the self-consistency results presented in the 
Appendix \^ we will write the effective phase action in 
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the following form 



S C fl#] = S a + (S) s { aa - bb \ 



( S2 )s\a,a,b,b] ~ ( S )l\a,a,b,b 



(70) 



The first term in Eq. ([7T|) could be omitted in the calcu- 
lations, because it gives the constant contribution to the 
total phase action. The averages in the right hand side 
of Eq. (|70|) are given in the sense that is discussed in Ap- 
pendix [A] From the second term in Eq. (j?0")) only the first 
term in the Eq. (|30p survives. Because we are searching 
the excitonic pairing and phase coherence mechanisms 
between the electrons and the holes in the system, we 
concentrate on the third term in Eq. (|70[) proportional to 
the average of the second order action and especially on 
the terms, which are proportional to the mixed products 
of the hopping amplitudes of two solid state bands tt or 
tt. These terms are responsible for the phase stiffness. 
We can rewrite the action Eq. (j?0")) in the form 



1 



(71) 



where the action <SoM is given in the Eq. ([31]) . The de- 
tailed calculation of the average of the second order ac- 
tion in Eq. (|7T|) is given in the Appendix [B] Keeping in 
evaluations only the terms proportional to the second 
order in the hopping amplitudes, we have for the phase 
stiffness action 



Sj M = -j 



dr cos 2 [np(rr) 

(rr'> 



-p(rV)], 



(72) 



where J is the phase stiffness parameter. The evalua- 
tions of J are given in details in the Appendix [B] and 
here we present only the phase stiffness critical lines. On 
the Fig. [5] we presented the phase stiffness parameter J 
for two different values of the hopping amplitude t. As 
the figures show, the phase stiffnesses (see the Appendix 
|B|) are strictly positive for all values of the Coulomb 
interaction. The upper curve in Fig. |6] corresponds to 
t = —0.3 and for the lower curve we have chosen the value 
t = —0.1. The phase stiffness is negligibly small, but is 
present along the whole EI region and it is especially 
important for the intermediate values of the Coulomb 
interaction, where the system of excitons undergoes the 
Bose-superfluid condensation. 



B. The condensate order 

Following discussion provided in the Section IIV1 the 
excitonic gap parameter at the very low temperatures 
(when the on-condensate excitations are negligible) could 
be written as in Eq. (1501) with the help of the condensate 



i 

O 




U/t 



FIG. 6: The zero temperature phase stiffness parameter J 
given in Eq. (|72[) for the three-dimensional system as a func- 
tion of the dimensionless interaction parameter U/t and for 
two different values of the hopping parameter t. 



transition probability amplitude ipo. In this section we 
will derive the excitonic condensate transition tempera- 
ture in the frames of the quantum rotor model. In the 
discussion above, we derived the form of the effective fr- 
action <S c ff [<p] = Sq [<p] + Sj [ip] . Here we start with the 
partition function given by 



Z = J [Dip] e 



(73) 



To proceed, we replace the phase degrees of freedom by 
the complex unimodular field z(rr) 



z{yt) = e 



(74) 



The variables z(rr) satisfy also the periodic boundary 
conditions z(r/3) — z(r0). This is done by implement- 
ing these new complex variables with the Faddeev- Popov 
ghost-fields method^ and we can introduce into the par- 
tition function in Eq. (l73l) the following identity 



VzVzS 



E 



xS(z 



)5{z 




1. 



(75) 



We introduce then the spherical constraint for the uni- 
modular complex variables z(yt) (|.z(rT)| 2 = 1). Indeed, 
the Eq. |75)) implies the evident condition 



N. 



(76) 
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The constraint of unimodularity can be resolved by intro- 
ducing the Lagrange multiplier A and using the functional 
representation of the Dirac delta-function. This adds a 
quadratic term (in the z-field) to the phase action 



^>(rr)| 2 -A 



xe 



-i/<f dT£„A( \z(r T )\ 2 ~l ) 



(77) 



Such a description is justified in the sense that the pair 
condensate has the macroscopic quantum phase coher- 
ence with the order parameter i/jq (see the section IIVI) . 
The phase stiffness action in Eq. (|T2"|) can be rewritten in 
more convenient form using the trigonometric transfor- 
mation formula 

cos 2 [<p(rr) — <p(t't)] = 2 cos 2 [p{rr) — (^(r'r)] — 1. 

(78) 

Then, in terms of the complex variables z(rr), the trans- 
formation in Eq. (|T5]) leads to a biquadratic form of the 
phase stiffness action. We have 



>j [z, z\ = 



J 
2 



/ dr ^ \z(xt)z{v't) + c.c.]' 



(79) 



This biquadratic action could be transformed into a 
quadratic one with the help of the Hubbard-Stratanovich 
transformation for the one component bose field, by in- 
troducing the complex variables w(yty't) at each pair of 
lattice positions r, r' and time r 

e i Si dT^ {rr , } [z(rr)z(r'r)+c.c.} 2 

JO aT ^{rr>) 2.1 x 



(80) 



= / [Vw]e 

S(f dTj2(rr>) W{TTT'T) [z (tT ) Z (V t) + C . C .] 



xe 



The saddle-point evaluation of the obtained action gives 
the critical value for the decoupling field w(rrr'r) 



w — J ([z(tt)z(t't) + c.c.}) = 2Jg, 



(81) 



where we used the expression of the bandwidth renor- 
malization factor g given in the Eq. (j42|) . Putting Eq.([8Tj) 
back into the expression Eq. ([5Dl we will have the phase 
stiffness contribution to the partition function in the 
spherical model 



4sJ£(„'> Si dT-z(vr)z(r'r) 



(82) 



Hence, the effective phase stiffness action Sj[z,z], after 
the decoupling procedure, is reduced to a quadratic one 
in z(yt) field, i.e. 



Sj[z,z]=-igjY^ / drz(vT)z{r' 

(rr>) J ° 



(83) 



Integrating the phase variables <p in Eq. ([73|) . we will have 
for the partition function of our quantum rotor approach 
the following expression 





p-\-ioo 


~v\~ 
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X 






J —ioo 


_2iri_ 




where Zq 



Z = Z J [DX] J [VzVz 



,— S\[z,z\ 



Z Q = J [Dip] e 



-Solv] 



(84) 



(85) 



S\ [z, z] in Eq. (|84|) is the total spherical action for our 
model. The form of the action S\ [z, z] in the Fourier 
space is 



S\[z, z] = — ^2 z(kuj n )G z 1 {kuj n )z(kuj n ) 



where 



f3N 
-PN\, 



g-\ku n ) = A - 8<?Je(k) + 7 _1 K) 



(86) 
(87) 



is the inverse of the bosonic charge propagator and u> n in 
Eq. ([86|) are the Bose-Matsubara frequencies. 7~ 1 (w n ) is 
the inverse of the Fourier transformation of the two-point 
phase correlator 7(rr, r'r') appearing after the phase in- 
tegration. It is defined as 

7 (rr,rV') = — / [Dip] e ~s [ v ]^[ v (rr)- v ( r 'r')] (88) 
Zo J 

The calculation of the Fourier transformation 7(w n ) of 
the function in Eq.(|88p is straightforward 2 ^ 



UZ Q 



E 



m— — oo 



1-16 [fy- 



where 



p \ 2 

u 1 



(89) 



(90) 



and the summation is over the winding numbers m, which 
define a manifold composed of an infinite set of the topo- 
logically equivalent paths. Along the critical line the in- 
verse of the uniform static order parameter susceptibility 
vanishes 



e z - 1 (k = o, w „ = o) = o, 



(91) 



which fixes the value of the Lagrange multiplier to the 
value A = Aq given by the equation 



A - 8.gJe(k = 0) + 7~V« = 0) = 0. 



(92) 



At the condensate phase transition boundary, we have 
the equation for the critical temperature 



— Y 



1 



PN^ A -8 5 Je(k)+7-i(O 



1. 



(93) 
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Considering the principal m = topological path, which 
is essential for the Bose-superfluid condensation, and 
performing the Bose-Matsubara frequency summation in 
Eq . (p?3"|) . we will have the following self-consistent equa- 
tion for the critical temperature 



-E 

2N ^ 



np (Cik) ~ n B (Cik) 
V/2 2 + 16[/gJ[e(0)-e(k)] 



= 1, 



(94) 



where tib (e) is the Bose-Einstein distribution function 
ns (e) = 1/ — l) and the variables Cik and ^2k are 
introduced as 



C2k = - 



-p + yjp 2 + WUJg [e(0) - e(k)] 



-p. - y/p? + IWJg [e(0) - e(k)] 



(95) 
(96) 



e(k) stands for the lattice dispersion for the three- 
dimensional cubic system and is given in Eq. (|4ip . The si- 
multaneous solution of the equation Eq. with the sys- 
tem of self-consistent equations Eq. (|6"Tj) from the Section 
IIVI gives the temperature transition critical line to the 
excitonic condensate state (EC) for the Falicov-Kimball 
model. The result is presented on the Fig. [7] (the lower 
curve on the figure). The upper curve in the panel rep- 
resents the temperature transition to the excitonic pair 
state (EP) for t = —0.3. This transition line defines the 
phase boundary at which the excitonic pair formation 
starts. Above the transition line the system is in the 
normal Fermi Liquid 25 (FL) state (N). Passing from the 
strong to the weak interaction regime, we observe that 
the macroscopic condensate begins its formation practi- 
cally at the same time with the excitonic pair formation 
and vanishes in a very weak interaction region (nearly 
at the same time that the exciton gap order disappears). 
For t = —0.1 we have presented the solution on the Fig. [5] 
We see that with decreasing the hopping parameter for 
the initial /-orbital electrons we obtain the similar tem- 
perature transition lines but with lower transition am- 
plitudes, so diminishing the hopping energy has a de- 
structive influence on the condensate transition and for 
a very low hopping amplitudes we could have the total 
destruction of the condensate state. From this point of 
observations, we can suggest that the condensate phase is 
due to the itinerant electronic correlations and localizing 
them leads to the destruction of the excitonic condensate 
state. It is clear that even very small phase correlations 
as exist in our model lead to the comparatively large (in 
amplitude) phase transition into the condensate regime 
with non-zero critical temperature. The maximum of 
the condensate transition coincides with the maximum 
of the phase stiffness plot (see Fig. [5]) and this is due to 
the fact that both transitions are strongly related. From 
the behavior of the excitonic gap parameter (Fig. [2]) it 
is clear that the transition maximum is displaced from 
that of the phase stiffness (Fig. [5]) and the condensate 
order (Fig. [7]). The maximum for the gap (in the case of 




FIG. 7: The transition temperature lines in the three- 
dimensional Falicov-Kimball model for t = —0.3 as a function 
of the normalized Coulomb interaction strength U. The upper 
curve represents the excitonic pair (EP) transition tempera- 
ture line Ta i.e. when A = and i^o = 0. The lower curve 
represents the excitonic condensate (EC) transition tempera- 
ture T c with A / and i/) / 0. 



t = —0.3) is achieved at U/t ~ 7.8, which corresponds to 
an imbalance of the c and / band electron concentrations 
nb — n a = 0.8 (this is near to the insulator regime), while 
the maximums for J and T c are pointed nearly at the 
same value of U/t = 5.1, which corresponds to the limit 
of semiconductor with nearly equal values of the electron 
concentrations n a rj nb- Physically, this is due also to 
the fact that the excitonic pair formation begins at the 
same time that the phase coherence appears in the sys- 
tem but the maximum number of excitonic pairs in the 
system is reached in the strong interaction region, while 
for the condensed pairs the maximum is shifted to the 
small and medium interaction regions. Another interest- 
ing physical interpretation of this phenomenon is that 
the phase coherence in the system is responsible for the 
Bose-superfluid phase transition and might cause also the 
excitonic pair formation, but, in essence, they still uncor- 
related. 



VI. SUMMARY 

In the present paper we have shown the existence 
of the Bose-superfuild condensate order in the three- 
dimensional strongly correlated electron system. We con- 
sidered the phase stiffness and we derived the effective 
phase stiffness action in the model. On the base of the 
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FIG. 8: The transition temperature lines in the three- 
dimensional Falicov-Kimball model for t — —0.1 as a function 
of the normalized Coulomb interaction strength U. The upper 
curve represents the excitonic pair (EP) transition tempera- 
ture line Ta i.e. when A = and tpo = 0. The lower curve 
represents the excitonic condensate (EC) transition tempera- 
ture T c with A / and i/) / 0. 



phase coherence between the conduction band electrons 
and valence band holes, we derived the expression for 
the phase stiffness parameter and we have shown that it 
is positive along the whole pair formation region of the 
normalized Coulomb interaction parameter. We proved 
that the phase stiffness is important for the macroscopic 
phase coherence order in the system, leading to the con- 
densation of the preformed excitonic pairs into the fun- 
damental state of the reciprocal momentum state. From 
our results it is clear that the critical temperature of the 
excitonic pair formation is not coinciding with the Bose- 
superfluid critical temperature. The Bose condensation 
of the excitonic pairs appears at the temperatures much 
smaller than the excitonic pair formation critical temper- 
ature is. Another important result is that of the existence 
of the BCS-BEC type crossover. In the small interaction 
limit the excitonic pairing and the condensate critical 
temperature lines are strongly coinciding. This is consis- 
tent with the superconducting transition scenario of the 
BCS model, in which the Cooper-pair formation and its 
condensation occur simultaneously and the system is in 
the BSC-like regime. For the strong interaction case we 
have the separation of the condensate order and the ex- 
citonic pairing critical temperature lines and the system 
is formally in the BEC-like regime. As a possible eval- 
uation of our results in future would be the determina- 
tion of the single-particle excitation spectra and density 



of states in the frames of the EFKM model. This will 
give us an access to compare our theoretical results with 
the temperature-dependent photoemission spectroscopy 
measurements. 
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Appendix A: The method of effective actions 

1. Fermionic action 

We start here with the partition function given in the 
Eq . (f3"2)l . We treat our problem of strongly correlated 
electrons in a self-consistent-field approximation based on 
the U(l) gauge transformation introduced in the Section 
TTTJ This is in general equivalent to the second order in the 
Fermi- liquid theory^ First, we integrate out the phase 
variables (p in Eq.(| 



Z = / [VaVa] [DbVb] e ~ s ^ a ^ b \ 



(Al) 



where the effective fermionic action S[a, a, 6, b] in the ex- 
ponential is defined as 



S[a, a, b, b] = — In / [Dip] e 



-S[a.a,b,b,<p] 



(A2) 



Then we expand the logarithm keeping only the terms 
up to second order in the hopping amplitude. We get for 
the effective action 

S[a,a,b,b}= S + (S) S 



Here 



/[2ty], 



(A3) 



(A4) 



The phase averaging can be considered as the quantum 
statistical averaging with an effective phase field action 
S e s [if]. This consideration is sometimes called like the 
self-consistency condition. Thus, we replace 



and 



_ J [Dip] ...e- 5rff M 



Then, for the effective fermionic action we have 



(A5) 



(A6) 



S[a,a,b,b]= S + <<5> SrffM 



(A7) 
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2. Phase action 

The integration over the fermions gives the effective 
action for the phase sector discussed above. The partition 
function in this case is 



Z = J [Dcp] , 



(A8) 



where the effective phase action in the exponential is 



S e s[<p] = - In / [DaVa] [DbVb] e" 



(A9) 



Again, expanding the logarithm we will have 

SoffM = s a + (s) ab 



1 



Here the fcrmionic average {■■■} ab is 
/ [DaVa] [VbVb] 



,-s 



(■■■) ab 



J [DaVa] [VbVb] e~ s 



(A10) 



(All) 



Repeating the same steps as in the previous case and 
using the self-consistency condition 

(S) ab -> {S) s[n , a ,l, b ] (A12) 
we get for the effective phase action in this case 



( S2 ) S\a,a,b,b] ~ ( S )s\a,a,b,< 



(A13) 



where Sq is an unimportant constant. After the self- 
consistency assumption, the average in Eq. (|Alll) will be 



')<S[a,a 



b.b 



J [DaVa] [VbVb] _ e s[a,a,b,b] 
/ [VaVa] [VbVb] e 



-S\a,a,b.b\ 



(A14) 



The relations in Eq. (IA6l) and Eq. (IA14p are princi- 
pal for the theory. Using these relations together, 
we resolve the problem of interacting fermions as a 
self-consistently coupled problem of the noninteracting 
fermions and the gauge-bosonic phase field. The bosonic 
sector plays the role of a suitable background, on which 
the collective excitations in the fermionic sector appear. 



Appendix B: Phase stiffness 

Here, we present a fermion averaging routine for av- 
eraging the four-fermionic type terms appearing in the 



expression of the third term of the effective phase ac- 
tion given in Eq. tfTTj) . This method is called the Wick 
averaging procedure and is used in many quantum-field 
theoretical treatments (see for example in Ref. |24| ) . We 
present here in details the derivation of the term propor- 
tional to tt. The other term proportional to tt can be 
derived analogically. For the fermion orbitals the mixed 
term proportional to tt is 

-\ E E I dTdT ' W^'i)t{^v' 2 )x 

(riri) {r2r' 2 ) 

<a(r 1 T)a(r' 1 r)fe(r 2 T')6(r^r')) x 

xe — i [v > (riT)-v(.r' 1 T)] e -i[ V (r 2 T')-<j>(r' 2 T')] 

+t(r 1 r' 1 )t(r' 2 r 2 ) <a(r 1 T)a(r' 1 r)6(r^')&(r 2 r')) x 

xe - i [v>{riT)-<p(r' 1 T)] e i[(p(r 2 T')- v >(r' 2 T')] 

+i(rin)t(r 2 r' 3 ) (a(rir)a(riT)6(r 2 T / )6(r / a T / )) x 

xe i[ V (r 1 T)-<p(r' 1 T)] e -i[<p(r 2 T , )-ip(r' 2 T')] 

+t(r' 1 r 1 )t(r' 2 r 2 ) (aftT)a(r 1 T)6(r / a T>(r a T')) x 

xe i[(f(T 1 T)-ip(r' 1 T)] e i[tp(r 2 T')- V >(r' 2 T')] ^ 



As an example, we give the Wick averaging result of the 
first four- fermion term in the expression of Eq. (|Bip 



<a(r 1 r) a (r' 1 r)6(r 2 r')6(r 2 r')) = 
(a(nr) a (rir)) {b(r 2 r')b(r' 2 r')) 

(a(r 1 r)6(r 2 r'))(a(r' 1 r)6(r 2 r')) 

= g aa (r' 1 -r 1 ,0)g bb (r 2 -r 2 ,0) 

-g ab (n - r' 2) r - T>)G ba {r 2 - r' 1( r' - r). (B2) 

We will keep only the terms proportional to the elec- 
tron density type order and the excitonic gap order. We 
neglect the other term containing the superconducting 
type ordering (a(rT)b(r' t')) or (a(rr)&(r'r')) because of 
the symmetry of the action in Eq. (|36[) . In turn, for the 
first treatment, the terms containing the density type or- 
derings (a(rr)a(r'r')) and (6(rr)6(r'r')) could be also 
omitted, since they are not contributing to the excitonic 
pair formation and condensate order. After calculating 
all averages in Eq. (IBl[) and recombining them with the 
similar terms coming from the component proportional 
to tt, we obtain the phase stiffness action in the form 
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Sj 



I dr {J(rr, r'r') cos [<p(tt) + <^(rr') — ^(r'r) — (p(r'r')] 
Jo 

+G ab (0, t - T')g ba (0, t' - t) cos [(p(rr) - p(r/) - <^(rV) + ^(rV)]} 

I 



(B3) 



where we introduced the interaction parameter J as 

J = -4tt / dr'g ab (r - r',r - r')e 6a (r - r',r' - r). 
Jo 



(B4) 



much faster than in the arguments of the trigonometric 
functions multiplied to it. For this reason, we will put 
t' = t in all arguments of the trigonometric functions 
corresponding to the zero order term in the Taylor series 
expansion of the phase variables 



</j(rr') = (p(rr) + (r' — r) d T (p(rr) + O (r' — r) 



The fermion propagators in Eq. (|B4[) are defined in the 
Section IIV1 The second term in Eq. (IB3|) gives the con- 
stant contribution to the phase stiffness action. The vari- 

(B5) 

ation of the product of the anomalous propagators with 

time difference |r — r' in the first term in Eq. (|B3[) is For the product of the anomalous propagators we have 
I 



Gab(r - r')G ha (r' - r) 



Z 2 ((3N) 2 LsL^ I AP] [SZMK'M) - |A| 21 



EE 



e(k)e(k')e- 



i(y n -w n )i5 



(B6) 



r 



Here z = 6 is the number of the nearest neighbor sites subara summations over the fermionic Matsubara fre- 
on the three-dimensional cubic lattice. After integrating quencies and we will have from the Eq. (IB4p and Eq. (IB6| 
over the imaginary time in Eq. (lB4[) . we perform the Mat- the phase stiffness parameter 



_ 16A 2 tt ^ e (k) e (k') 



2 2 A^ 2 



kk' 



Ai(kk') tanh ( - A 2 (kk') tanh 



r 



(B7) 



The parameters Ai(kk') and A2(kk') entering in Eq. (|B7l 
are given by 



Ai(kk') = 
A 2 (kk') = 



1 



E+(k) - E+(k>) E+QK)-E-Qs?y 
1 1 



E-(k)-E-(k') E-(k)-E+(k') 



(B8) 



Following the Section IIV1 the summations over the k- 
wave vectors in Eq. (|B7J) could be transformed into the 
integrations with the help of the appropriate density of 
states given in Eq. (l67[) . The numerical evaluation results 
of the expression Eq. (IB7p for T = K are presented in 
the Fig. |6]and discussed in the Section [V] of the present 
paper. 
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